Effects of hadronic potentials on elliptic flows in relativistic heavy ion collisions 
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Within the framework of a multiphase transport (AMPT) model that includes both initial partonic 
and final hadronic interactions, we show that including mean-field potentials in the hadronic phase 
leads to a splitting of the elliptic flows of particles and their antiparticles, providing thus a plausible 
explanation of the different elliptic flows between p and p, K + and K~ , and 7r + and ir~ observed 
in recent Beam Energy Scan (BES) program at the Relativistic Heavy-Ion Collider (RHIC). 
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Understanding the phase diagram of strongly inter- 
acting matters is one of the main goals of the experi- 
ments that are being carried out in heavy-ion collisions 
at RHIC For collisions at its highest center-of-mass 
energy of y^sivjv = 200 GeV, convincing evidences have 
been established that a strongly interacting quark-gluon 
plasma (sQGP) is formed in these collisions 0. The 
baryon chemical potential of the formed sQGP is small, 
and according to calculations based on the lattice quan- 
tum chromo dynamics (QCD), the phase transition be- 
tween the sQGP and the hadronic matter for such a small 
baryon chemical potential is a smooth crossover transi- 
tion Q . Theoretical studies based on various models have 
shown, on the other hand, that this crossover transition 
would change to a first-order phase transition when the 
baryon chemical potential is increased [l|. To search for 
the critical point in the baryon chemical potential and 
temperature plane at which such a transition occurs, the 
BES program involving Au+Au collisions at lower ener- 
gies of ^snn = 7.7, 11.5, and 39 GeV has recently been 
carried out. Although no definitive conclusions on the 
location of the critical point in the QCD phase diagram 
have been obtained from these experiments, several inter- 
esting phenomena have been observed Q ■ Among them 
is the smaller elliptic flows of p, K~ and 7r+ than those 
of p, K + and ir~ , respectively. The differences decrease 
with increasing collision energy and become essentially 
a constant at higher collision energies These sur- 
prising results were recently attributed to the different 
elliptic flows of transported and produced partons dur- 
ing the initial stage of heavy ion collisions @. Also, it 
was suggested that the chiral magnetic effect induced by 
the strong magnetic field in non-central collisions could 
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be responsible for the observed difference between the 
elliptic flows of ir + and n~ Q. 

It is known from heavy ion collisions at lower colli- 
sion energies at SIS/GSI and AGS/BNL that the elliptic 
flow of nucleons is affected not only by their scattering 
but also by their mean-field potentials in the hadronic 
matter Q. In particular, in collisions at a few AGcV 
when participating nucleons are not blocked by specta- 
tor nucleons, an attractive (repulsive) potential is found 
to result in a smaller (larger) elliptic flow. This is because 
particles with attractive potentials are more likely to be 
trapped in the system and move in the direction perpen- 
dicular to the participant plane while those with repulsive 
potentials are more likely to leave the system and move 
along the participant plane, thus reducing and enhanc- 
ing their respective elliptic flows. Also, the potentials of 
a particle and its antiparticle are different, and they gen- 
erally have opposite signs at high densities [j| [Ic} ■ As 
a result, particles and antiparticles are expected to have 
different elliptic flows in heavy ion collisions when the 
produced matter has a nonzero baryon chemical poten- 
tial. Furthermore, the difference between the potentials 
of a particle and its antiparticle diminishes with decreas- 
ing baryon chemical potential, so their elliptic flows are 
expected to become similar in higher energy collisions 
when more antiparticles are produced. These effects are 
all consistent with what were seen in the experimental 
data from the BES program. 

In the present paper, we study quantitatively the el- 
liptic flows of particles and their antiparticles at BES en- 
ergies by extending the AMPT model [ll| to include the 
potentials of baryons, kaons, and pions as well as their an- 
tiparticles. The AMPT model is a hybrid model with the 
initial particle distributions generated by the Heavy-Ion 
Jet Interaction Generator (HIJING) model [l2| via the 
Lund string fragmentation model. In the string melting 
version of the AMPT model, which is used in the present 
study, hadrons produced from excited strings in the HI- 
JING model are converted to their valence quarks and 
antiquarks, and their evolution in time and space is then 
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modeled by Zhang's parton cascade (ZPC) model [llj]. 
Different from previous applications of the AMPT model 
for heavy ion collisions at higher energies, the parton 
scattering cross section and the ending time of the par- 
tonic stage are adjusted in the present study to ap- 
proximately reproduce measured elliptic flows and the 
hadronic energy density (~ 0.30 — 0.35 GeV/fm 3 ) at the 
extracted baryon chemical potential and temperature at 
chemical freeze out [ijj]. Specifically, we take the parton 
scattering cross section to be isotropic with the value 3, 6, 
and 10 mb and the ending time of the partonic stage to be 
3.5, 2.6, and 2.9 fm/c for collisions at y%iv = 7.7, 11.5, 
and 39 GeV, respectively. At hadronization, quarks and 
antiquarks in the AMPT model are converted to hadrons 
via a spatial coalescence model, and the scatterings be- 
tween hadrons in the hadronic sta ge a re described by a 
relativistic transport (ART) model [15| that has been ex- 
tended to also include particle-antiparticle annihilations 
and their reverse reactions. 

For the nuclcon and antinuclcon potentials, we 
take them from the relativistic mean-field model used 
in the Relativistic Vlasov-Uehling-Uhlenbeck transport 
model [3, that is 

U n ,n(pb,P§) = ^s{pb,Pb) ± ^Z{pb,Pb)^ (!) 

in terms of the nuclcon scalar E s (ps, Pb) an d vector 
Y*®(pb,Pb) self-energies in a hadronic matter of baryon 
density ps and antibaryon density p§. The "+" and "— " 
signs are for nucleons and antinuclcons, respectively. We 
note that nucleons and antinucleons contribute both pos- 
itively to £ s but positively and negatively to £„, respec- 
tively, as a result of the G-parity invariance. Since only 
the light quarks in baryons and antibaryons contribute to 
the scalar and vector self-energies in the mean-field ap- 
proach, the potentials of strange baryons and antibaryons 
are reduced relative to those of nucleons and antinucleons 
according to the ratios of their light quark numbers. 

The kaon and antikaon potentials in the nuclear 
medium are also taken from Ref. [l6| based on the chiral 
effective Lagrangian, that is U K g = uj k g — uiq with 

uigg = sjm\ +p 2 - aggp s + (6jrf) 2 ± b K ff§{2) 

and ujq = y m /<- + P 2 i where vtlk is the kaon mass and 
a K = 0.22 GeV 2 fm 3 , a R = 0.45 GcV 2 fm 3 and b K = 0.33 
GeVfm 3 are empirical parameters taken from Ref. fl7j |. 
In the above, p s is the scalar density, which can be deter- 
mined from ps and pg through the effective interaction 
used for describing the properties of nuclear matter, and 
PB Ct = Pb — Pb 1S the net baryon density. The " +" and 
"— " signs are for kaons and antikaons, respectively. 
The pion potentials are related to their self-energies 
according to U^±o = n^°/(2m 7r ), where is the 
pion mass. In Ref. [la ], the contribution of the pion- 
nucleon s-wave interaction to the pion self-energy has 
been calculated up to the two-loop order in chiral pertur- 
bation theory. In asymmetric nuclear matter of proton 



density p p and neutron density p n , the resulting n and 
7r + self-energies are given, respectively, by 

nj(p P ,p„) = p n [T~ N -T+ N ] -p P [T~ N +T+ N ] 

+ n i ~ol(ft" Pn) + n cor(Pp> Pn) 

U+(p p , Pn ) = U-(p n ,p p ). (3) 

In the above, are the isospin-even and isospin-odd 
ttN s-wave scattering T-matrices, which are given by 
the one-loop contribution in chiral perturbation theory 
and have the empirical values T~ N 1.847 fm and 
^ttN ~ — 0.045 fm extracted from the energy shift and 
width of the Is level in pionic hydrogen atom; ±I~[ is 
due to the relativistic correction; and n^ r is the con- 
tribution from the two-loop order in chiral perturbation 
theory. Their expressions can be found in Ref. [HI . For 
nucleon resonances and strange baryons in a hadronic 
matter, we simply extend the above result by treating 
them as neutron- or proton-like according to their isospin 
structure [l5T | and light quark numbers. Because of the 
G-parity invariance, the contributions of antiprotons and 
antineutrons in the hadronic matter are similar to those 
of neutrons and protons, respectively. We neglect in the 
present study the pion- nucleon p-wave interaction [l9| . 
which is expected to reduce and enhance, respectively, 
the 7r + and ir~ potential difference due to the s-wave in- 
teraction in the neutron-rich and proton- rich matter [20| , 
since its inclusion in the transport model is highly non- 
trivial (u. 




FIG. 1: (Color online) Mean-field potentials of N (a), K (c), 
and K (d) at zero momentum as functions of baryon and 
antibaryon densities and of n~ (b) as a function of neutron- 
like and proton-like densities. 

In Fig. [TJ we show the N, K, and K potentials as 
functions of ps and p§ and the ir~ potential as a func- 
tion of neutron- like and proton- like densities p n and p p . 
The N and ir + potentials are related to those of TV and 
tt~ by Un(p b ,Pb) = U N (p B ,p B ) and U„+{p p ,p n ) = 
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U 7T -(p n , p p ). In the absence of antibaryons, the N po- 
tential is slightly attractive while that of N is strongly 
attractive, with values of about —60 MeV and —260 
MeV, respectively, at normal nuclear matter density 
po = 0.16 fm~ 3 . The latter is similar to that deter- 
mined from the non-linear derivative model for small 
antinucleon kinetic energies [22[ and is also consistent 
with those extrapolated from experimental data [23T - i28| . 
For pions in neutron-rich nuclear matter, the potential 
is weakly repulsive and attractive for 7r~ and 7r + , re- 
spectively, and the strength at po and isospin asymmetry 
S = (p n — Pp)/(Pn+Pp) = 0.2 is about 14 MeV for 7r _ and 

— 1 MeV for ir + . In antibaryon-free matter, the K poten- 
tial is slightly repulsive while the K potential is deeply 
attractive, and their values at po are about 20 MeV and 

— 120 MeV, respectively, similar to those extracted from 
the experimental data [29|-[32| and used in the previous 
study [33l |. 
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FIG. 2: (Color online) Time evolution of the baryon den- 
sity (a), antibaryon density (b), neutron- like density (c), and 
isospin asymmetry (d) in the central region of the hadronic 
phase from the AMPT model for Au+Au collisions at b — 8 
fin and v / sjvjv = 7.7, 11.5, and 39 GeV. 

Figure [2] displays the time evolution of the baryon den- 
sity (a), antibaryon density (b), neutron-like density (c), 
and isospin asymmetry (d) in the central region of the 
hadronic phase in Au+Au collisions at impact param- 
eter b = 8 fm for the three different BES energies of 
y/sNN = 7.7, 11.5, and 39 GeV. It is seen that the baryon 
density decreases while the antibaryon density increases 
with increasing collision energy, resulting in a decrease of 
the net baryon density with increasing collision energy. 
Also, the isospin asymmetry is very small in the hadronic 
phase for all three energies due to the considerable num- 
ber of A hyperons which do not carry isospin and the 
larger number of n~ than tt + produced in the collisions. 

The differential elliptic flows of p, K + , and 7r + as well 
as their antiparticles with and without hadronic poten- 
tials at three different BES energies from the string melt- 
ing AMPT model are shown in Fig. [31 They are calcu- 
lated with respect to the participant plane, that is Vi = 
(cos[2(r/>— ^2)]), where <fi = atan2(p y , p x ) is the azimuthal 
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FIG. 3: (Color online) Differential elliptic flows of mid- 
rapidity (\y\ < 1) p and p [(a), (b), (c)], K+ and K~ [(d), 
(e), (f)], and 7r + and tt~ [(g), (h), (i)] with and without 
hadronic potentials U in Au+Au collisions at b = 8 fm and 
y'sjvat = 7.7, 11.5, and 39 GeV from the string melting 
AMPT model. 



angle, * 2 = [atan2((rp sin 2r/> p ), (r£ cos 2<j> v )) +tt]/2 is the 
angle of the participant plane, and r p and <f> p are the 
polar coordinates of the participants. Without hadronic 
potentials the elliptic flows from the AMPT model are 
similar for particles and their antiparticles. Including 
hadronic potentials increases slightly the p and p elliptic 
flows at pt < 0.5 GeV/c, while reduces slightly (strongly) 
the p (p) elliptic flow at higher pr, consistent with the 
expectations from the relative strength of the attractive 
potentials for N and N shown in Fig. [TJ Hadronic poten- 
tials also increase slightly the elliptic flow of K + while 
reduces mostly that of K~ , again consistent with what is 
expected from the K and K potentials shown in Fig. [I] In 
addition, the effect from the potentials on the elliptic flow 
decreases with increasing collision energy, which is consis- 
tent with the decreasing baryon density and net baryon 
density of produced hadronic matter with increasing col- 
lision energy shown in Fig. [5] The difference between the 
differential elliptic flows of p and p, and between those of 
K + and K~ below y / 5jv r /v = 11.5 GeV are qualitatively 
consistent with the experimental data Q, while that of 
7r~ and 7r + is small in all three energies due to the small 
isospin asymmetries shown in Fig. [21 

Our results for the relative pr-integrated «2 differ- 
ence between particles and their antiparticles, defined by 
[v2(P) — V2(P)]/v2(P), with and without hadronic po- 
tentials arc shown in Fig. 2) These differences are very 
small in the absence of hadronic potentials. Including 
hadronic potentials increases the relative V2 difference 
between p and p and between K + and K~ up to about 
30% at 7.7 GeV and 20% at 11.5 GeV but negligibly at 
39 GeV. These results are qualitatively consistent with 
the measured values of about 63% and 13% at 7.7 GeV, 
44% and 3% at 11.5 GeV, and 12% and 1% for the rela- 
tive V2 difference between p and p and between K + and 
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FIG. 4: (Color online) Relative elliptic flow difference between 
p and p, K + and K~ , and ir + and it" with and without 
hadronic potentials U at three different BES energies from 
the string melting AMPT model. Results for different species 
are slightly shifted in energy to facilitate the presentation. 



K~ , respectively Similar to the experimental data, 
the relative v 2 difference between n + and 7r~ is negative 
at all energies after including their potentials, although 
ours have smaller magnitudes. We have also found that, 
as seen in the experiments Q, the relative v 2 difference 
between A hyperons and A is smaller than that between 
p and p, because the A (A) potential is only 2/3 of the 
p(p) potential. 

To summarize, we have studied the elliptic flows of 
p, A" + , 7r + and their antiparticles in heavy ion colli- 
sions at BES energies by extending the string melting 
AMPT model to include their mean-field potentials in 
the hadronic stage. Because of the more attractive p than 
p potentials, the attractive K~ and repulsive K + poten- 
tials, and the slightly attractive 7r + and repulsive 7r~ po- 
tentials in the baryon- and neutron-rich matter formed 
in these collisions, smaller elliptic flows are obtained for 
p, K~, and ir + than for p, K + , and tt~ . Also, the dif- 
ference between the elliptic flows of particles and their 



antiparticlc is found to decrease with increasing collision 
energy as a result of decreasing baryon chemical potential 
of the hadronic matter. Although our results are qual- 
itatively consistent with the experimental observations, 
they somewhat underestimated the relative elliptic flow 
difference between p and p as well as that between n~ and 
7r + and overestimated that between K + and K~ . In our 
studies, we have, however, not included other effects that 
may affect the v 2 difference between particles and their 
antiparticles. For example, we may have overestimated 
the annihilation between baryons and antibaryons as this 
could be screened by other particles in the hadronic mat- 
ter HH . Including the screening effect would increase the 
duration of the attractive potential acting on antibaryons 
and thus reduces their elliptic flow, leading therefore to 
an increase in the difference between the elliptic flows 
of baryons and antibaryons. Also, the different elliptic 
flows between particles and their antiparticles are as- 
sumed in the present study to come entirely from the 
hadronic mean- field potentials. As shown in Ref. (35j . 
the collective flow of partons can also be affected by their 
mean-field potentials in the partonic matter. If quarks 
and antiquarks have different mean-field potentials in the 
partonic matter, this would then lead to different ellip- 
tic flows for particles and their antiparticles in the initial 
stage of the hadronic phase after hadronization. It will 
be of great interest to include in future studies these ef- 
fects as well as the effect due to different elliptic flows 
between produced and transported partons [6j and the 
chiral magnetic effect Q in order to understand more 
quantitatively the different elliptic flows between parti- 
cles and their antiparticles observed in relativistic heavy 
ion collisions. 
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